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We examine probability distributions for thermodynamic quantities in finite-sized random sys- 
tems close to criticality. Guided by available exact results, a general ansatz is proposed for replicated 
free energies, which leads to scaling forms for cumulants of various macroscopic observables. For the 
specific example of a planar flux line lattice in a two dimensional superconducting film near Hci , we 
provide detailed results for the statistics of the magnetic flux density, susceptibility, heat capacity, 
and their cross-correlations. 
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Impurities in a sample are expected to modify vari- 
ous measurements, making it desirable to characterize 
the probability distribution functions (PDFs) for the out- 
comes. These PDFs may provide important insight about 
the underlying physics; as in the case of universal con- 
ductance fluctuations in mesoscopic circuits, a subject of 
much recent investigation Q|. Here we consider the sig- 
nature of impurities in thermodynamic systems at equi- 
librium. At one extreme, microscopic quantities such as 
two-point correlation functions are quite sensitive to dis- 
order; in some cases their PDFs exhibiting complicated 
multiscaling behavior On the other hand, the free en- 
ergy, and other macroscopic properties, are expected to 
be self-averaging, converging to fixed values in an infinite 
system. If all correlation lengths in the system are finite, 
the PDFs for a mesoscopic (finite-sized) sample should 
be governed by the central limit theorem. We thus focus 
on systems with long-range correlations, such as close to 
a critical point, or in a flux line (FL) lattice. 

In the most interesting cases, disorder is relevant, lead- 
ing to novel correlations distinct from the pure case. 
However, due the difficulty of characterizing collective 
behavior in such glassy disordered systems ^J, little is 
known about the corresponding PDFs. Aharony and 
Harris recently studied the PDFs of thermodynamic 
observables, near critical points with relevant random- 
ness, in finite-sized random systems. They find a lack of 
self-averaging, and universal, non-Gaussian PDFs. There 
are also some results for elastic manifolds pinned by im- 
purities: Using replicas, Mezard and Parisi relate the 
scaling of the PDF for susceptibility with size of a man- 
ifold, to its roughness exponent ||]. Exact results on 
PDFs have been obtained for a directed polymer on a 
disordered Cayley tree by Derrida and Spohn, using a 
mapping to a deterministic differential equation 

An important example of pinned elastic media is pro- 
vided by FLs in a superconductor with point impurities. 
Indeed, our study was motivated by the experiment of 
Bolle et al. on a 2d FL lattice oriented parallel to 
a thin micrometer-sized film of 2H-NbSe2 . Magnetic re- 



sponse measurements show interesting sample-dependent 
fine structure in B{H); a fingerprint of the underlying 
pinning landscape. Interestingly, the problem of 2d lines 
with impurities is amenable to an exact solution using 
replicas, which gives not only the quenched average but 
also cumulants of the free energy Q . We shall use this ex- 
act solution to motivate an ansatz for the scaling of PDFs 
in general disorder-dominated thermodynamic systems. 

In what follows, we first present a general ansatz for 
the scaling of the replicated free energy. The key assump- 
tion is treating the number of replicas as a scaling field. 
Consequences of this ansatz for the scaling of cumulants 
of thermodynamic observables in mesoscopic samples are 
then enumerated. We then re-examine the example of 
FLs in a 2d layer in more detail, proposing specific ex- 
perimental tests of the theoretical picture. 

The lack of translational symmetry in the presence of 
impurities is cured by calculating disorder averaged mo- 
ments, [Z'^], of the partition function. These moments 
then provide information about PDFs of the free energy, 
magnetization, or susceptibility. We thus focus on the 
scaling of F„ = — Tln[Z"], the free energy of n repli- 
cas of the system, with interactions resulting from the 
disorder average. In fact Fn can be determined exactly 
for a single elastic line, or a lattice of non-crossing lines, 
in 2d with point impurities 1^. While this is sufficient 
for the FL experiment [Q, we would like to propose a 
more general scaling ansatz for random systems. Con- 
sider np replicas of a system of size with p 
groups of n replicas subject to the same reduced tem- 
perature Tj and scaling field ipji for j — f The 
vector of scaling fields if> can consist of say external 
magnetic fields, or chemical potentials. Our scaling 
ansatz for the singular part of the free energy density, 

4*^) {n,T,iP) = (n6^ii , rb^/"^^ , t/^fo*/''" ) . (1) 

To include FL systems, we have allowed for anisotropic 
scaling, with a 'roughness exponent' C relating the (d— 1)- 
dimensional longitudinal (Z||), and the Id transversal (/_l) 
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scales, by l± ^ Z^. (For isotropic systems, we can set 
C = 1.) Compared to the standard scaling hypothesis for 
pure systems, the replica structure introduces n as an 
additional scaling field with dimension 6*11 = C0±- The 
exponent appears in the modified hyperscaling rela- 
tion 2 — a = l + C^ ^||), which follows readily from 
Eq. (|l|) by comparing the 0{n) terms. Besides its agree- 
ment with available exact results, a mean-field analysis 
of the 0^-model with Gaussian random fluctuations in Tc 
also supports our ansatz. 

Scaling of disorder averaged cumulants of the ther- 
modynamic observables can now be extracted from the 
above ansatz for ln[n^^j Z"(tj-, -0^)]. Here we enumerate 
the main results (details to appear elsewhere): 

(1) Free energy cumulants: For p = 1, we have 

,„[^.,,„| = f (ZIIEM.. ,2) 

where [. . .]c denotes the connected or cumulant disor- 
der average. Using a scaling factor proportional to the 
correlation length, b ^ ^ t^"^" , gives [F'''{t)]c ^ 
.j-i/|i(ti-i+c-pe||)^ Pqj. mesoscopic systems of size L± L^, 
close to criticality (^|| > L||) we must set 6 L||, result- 
ing in [FP]c ^ -/jjl^" . The exponent 6*11 thus characterizes 
all sample-to-sample fluctuations in the free energy. 

(2) Thermal averages, such as the magnetization or 
the number of 'particles', are first order derivatives of 
the free energy, which we will denote generally by {X) = 
—dF/dip\^=o. For a specific realization of disorder, the 
partition function can be expanded as 

Z(^) = Z(0)f;^^^(X^), (3) 

where (. . .) denotes the thermal average. Using this 
representation, one can show that the p^^ coefficient of 
the expansion of ln[Z"'(V')] with respect to nip is the 
connected average [{X)p]cT^p /pi. Applying the scaling 
hypothesis of Eq.(|l|), we obtain with 5 = 2 ~ a — (3 
the critical scaling behavior [(X)^']^ ~ t'^p, with /3p = 
pP -v\\{d-l + C){P - 1), and /3 defined by [{X)] - . 

(3) Response functions: Cumulants of the susceptibil- 
ityj X = {{^'^) ~ can be calculated in a similar 
way. We now consider p sets of n replicas, each with the 
same scaling field ipj, j = 1, . . . ,p. The coefficient of the 
term n'P(ijji ■ . . . ■ ipp)'^ of the expansion of the generating 
function ln[J^^^j^ Z"'{'ipj)] is given by the cumulant aver- 
age {2T)~P[xP]c- This gives a critical scaling of the form 
[X% - t"^", with 7p =p7 + i/||(d - f + C)(p- f). 

(4) Cross correlations, such as F^y — {XY) — {X){Y), 
where 1" is a derivative of the free energy with respect 
to another scaling field ip, are also of interest. For ex- 
ample, the cross correlation, Tem, of magnetization M 



and energy E, has been measured in numerical simula- 
tions, since it allows for a more accurate estimate of the 
exponent ratio a/v then a direct method The gen- 
erator of TxY is ln[nj=i ■^^"(V'ij V'j)]; with coefficients 

r^P[r^y]c corresponding to the terms n^ipii^i ■ ■ ■V'pV'p- 
Our scaling ansatz yields Tem ~ rP^I^-^^-'^nid-i+OiP-i) ^ 
There are also correlations between different susceptibil- 
ities. Along the lines presented above, it is possible to 
show that these cross-correlations can be generated from 
the expansion of ^n[Y[^^i Z"{i/jj)Z"{ipj)], and that the 
cumulant (2T)~^^'[(xx)^]c is the coefficient of the term 
n'^Ptpfipf . . .tppipp. Choosing for x and x the magnetic 
susceptibility and the heat capacity c, respectively, we 

get [(xc)P]c - T-27p-i'||(d-l+C). 

Finally, we compare with the results of Aharony and 
Harris on the relative cumulants Rp,x = [X^]c/[X]p, 
in a (/('^-model with random Tc- By a perturbative renor- 
malization group (RG), in which they assume a Gaussian 
distribution for randomness on all length scales, they ob- 
tain i?p_^ = p!2P~'^i?2,x^j foi' the magnetic susceptibility. 
Indeed, for all observables (including susceptibilities) we 
find Rp,x ^ ^2^' hut the exact coefficient cannot be 
obtained from the scaling ansatz. However, even for an 
originally Gaussian random Tc, higher then second cu- 
mulants are generated by RG transformations and yield 
additional (universal) contributions to the above coeffi- 
cient Assuming that randomness is relevant at the 
phase transition, we find for the second relative cumu- 
lant of any observable X that i?2,x ~ (6/L||)''~^(6'^/L±). 
Choosing a scaling factor b ^ ^|| if ^|| ^ and b ^ L|| 
in the critical regime ^|| 3> we reproduce the results 
ofRef. g. 

We now apply these results to the experimental study 
by BoUe et al. of a planar, randomly pinned vortex array 
penetrating a mesoscopic quasi 2-dimensional thin single 
crystal of 2H-NbSe2 with weak pinning 0. The crystal 
was glued onto a silicon micromachined mechanical res- 
onator. By measuring jumps in the resonant frequency 
caused by magnetic FLs entering the sample, the num- 
ber of lines was determined very accurately as a function 
of the applied magnetic field. Close to the lower criti- 
cal field Hci, the jumps were irreversible, indicative of 
the difficulty of FLs finding optimal pinning configura- 
tions. At higher fields, the increased line density should 
result in a smaller pinning length, enabling the FL lat- 
tice to find its optimal pinning state. Indeed reversible 
behavior is observed experimentally in this case. In both 
regimes, the response depends on the detailed configura- 
tion of pinning sites, the sample geometry, and the vortex 
interactions. Therefore, the observed discrete jumps in 
the constitutive relation, B(H), provide a fingerprint of 
the disorder in a specific sample. 

Previous theoretical and experimental work concen- 
trated mainly on the determination of the relation B(H) 
near Hd for different kinds of disorder [p!o[ . Fluctuations 
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in B{H) are also of interest, characterizing more clearly 
the non-trivial effects of randomness in this mesoscopic 
system. Therefore, we study entire PDFs for observables 
like magnetic susceptibility Xi heat capacity c, and 
correlations between them. First consider a single FL to 
introduce the basic notations. The transversal wander- 
ing of such a line is described by a trajectory x{y), which 
has an energy |0] 



n = 



dy 



(4) 



Here, g = (p^Hd/ {'^tt) is the elastic stiffness of the FL 
carrying a flux quantum (pQ. The random potential 
V{x,y) mimics point like pinning, and is characterized by 
the disorder average \V{x,y)V{x' ,y')\ = AS^{x — x')6{y — 
y'), with strength A and a short range function, of 
width of the in-plane coherence length ^. The dimensions 
of the sample are taken to be L along the field direction 
y, and W in the transverse direction x. 

A single FL can wander freely in the transverse di- 
rection to take advantage of the randomly distributed 
pinning centers. This leads to an anomalous growth of 
its displacement by Sx{L) ^ L^/^ In contrast, in a 
lattice of lines, the non-crossing condition is a strong re- 
striction for the possible configurations of each line. The 
results of Ref. § for such a lattice, only apply to FLs at 
temperatures T larger than T* = {gS^A^/^, the height 
of the smallest energy barrier due to the random poten- 
tial. Recently, Korshunov and Dotsenko generalized 
the results for a single line to the low temperature limit, 
T 4C T*, by using a replica interaction potential with a 
small but finite curvature, instead of the rectangular well 
which works in the high-T limit. The generalization of 
their solution to an array of lines is straightforward. The 
non-trivial part of the free energy of n replicas of the sys- 
tem, with a fixed number of TV = W/uq of lines at mean 
separation oq, can be summarized in both limits as B 



with the parameter = 1 for T > T*, and A: « T/T* 
for T T*. Here we have neglected the trivial contri- 
bution to the free energy which is linear in N, and Q is 
an analytic function |l4[ . 

To study the fluctuations in the number of FLs near 
Hci, we use a grand canonical description with a chemi- 
cal potential /i = ghL, where h = {H — Hci)/Hci is the 
reduced magnetic field. The PDF of N is characterized 
by the disorder averaged cumulants [NP]c in the number 
of FLs, which are given by times the coefficients of 
the terms {nfi)P in the expansion of ln[Z"(/i)], leading to 



kA 



2-p 



^(5-3p)/2^ 



(6) 



The moments of the magnetic flux density follow from 
B = (t)oN/{dW), where d is the thickness of the sam- 
ple. It is interesting to note that the second cumulant 
is universal, [N^]c ~ {W/L)h~^^'^, independent of tem- 
perature and disorder strength, in both the high and low 
temperature regimes. For T <^T*, all cumulants are in- 
dependent of T, but do depend on the disorder strength. 
In the thermodynamic limit WL — > oo, all cumulants 
[BP]c approach zero for p > 1. On the other hand, for 
mesoscopic systems, the divergence of [BP]c for ft, — > 
is stopped if the average distance oq between the FLs 
approaches the system size W. Deep in the glassy low 
density phase near Hd, we obtain from Eq. (j^ for the 
free energy cumulants 



T y 
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/j(5-p)/2 



(7) 



This agrees with our scaling ansatz with exponents C = 
2/3, 0|| = 1/3, J^ii — 3/2 and 6 — 1, and gives a universal 
value for p — 2, which still depends on the stiffness g. 

Now consider the response of the FL lattice to changes 
in the external magnetic field, measured by the suscep- 
tibility X- This is related to thermal fluctuations in the 
number of FLs by X = (L/M^)((iV2)-(iV)2)/r, for a fixed 
realization of disorder. The disorder averaged p'^ cumu- 
lant of the number fluctuations can be obtained from the 
generating function ln[J^^^-^ for p different chem- 

ical potentials iij , giving the susceptibility cumulants 



{WL)^-Pg^-'^P 



T 
kA 



2-p 



^-5(p-l)/2^ 



(8) 



This result deserves a few comments: First, the disorder 
averaged susceptibility [x] is non-singular at the transi- 
tion h —^ 0, since the exponent vanishes forp = 1. In fact, 
due to a statistical tilt symmetry ||l5| the susceptibility is 
simply related to the compression and tilt elastic moduli 
(cii, C44) of the FL lattice, by [%] — (27r/ao)^(ciiC44)^"^^^- 
This is in agreement with our result [x] T/{kA), as 
can be seen by using C44 — g/ao, rewriting cn in terms 
of the steric repulsion between the lines, and estimat- 
ing ao = W/N from Eq. (^) with p = 1, and a value 
for h of order one. Again, the variance is the only 
moment of the susceptibility, which shows universality, 
[X% ^ {WLg^r^h-^/^ !§. 

Next, consider the response of a fixed number of FLs to 
changes in temperature. Keeping the number N of FLs 
constant by adjusting the magnetic field, but changing 
the temperature, allows us to study thermal fluctuations 
around the global ground-state. Physically, the response 
of the FL lattice to changes in magnetic fleld H , or tem- 
perature T, are quite different: A change in H (or N) 
usually leads to a complete rearrangement of the whole 
ensemble of lines. Increasing T , however, causes stronger 
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entropic repulsions between the N lines, which now fluc- 
tuate around their state of optimal pinning. The response 
to a change in T also depends on the detailed pinning 
landscape, producing a heat capacity c which is sample 
specific. Therefore, the statistics of c are also of inter- 
est; its cumulants can be obtained from the generating 
function ln[]^^^-^ Z"{Pj)] with p different temperatures 
Tj = l//3j , as (for fixed N) 



p—3 P+1 J 

LW—{kA) — g- 



— N—T- 



(9) 



In the high temperature limit [T ^ T*), the moments 
[c^]c r-(3p+i)/2 decay faster with temperature than in 
the low temperature limit (T < T*), where [cP]c T-p. 

We may also examine cross-correlations between the 
different susceptibilities such as the heat capacity c, and 
the magnetic susceptibility x- In. the thermodynamic 
limit, i.e. beyond a characteristic system size, different 
susceptibilities are expected to be statistically indepen- 
dent. For example, this has been demonstrated for the 
magnetic susceptibilities of two noninteracting FL lat- 
tices with different random potentials The correla- 
tions between c and x can be obtained from our scaling 
theory by extracting the coefficients of . . . ^p/3p 

of ln[ni=i ^"(a^j) 11^=1 Using the results of 

Eqs. (|) and (|), we get 



[(Xc)P]e 



{wiy 



(10) 



For finite systems, the divergence as ft, is cut off at 
h ^ kA/ (gTW). Then the result can be rewritten as 
[(xc)^]c/([x^]c[c^]c) ^ Lc/L, with a characteristic length 
scale Lc — {gT/kAY^'^W'^^'^ for the decay of correlations. 
This length-scale has a simple physical interpretation: Lc 
is the length of a single FL whose transverse wander- 
ings can explore the whole sample width W. Since the 
transversal fluctuations of a line of length L is grow as 
Sx{L) « {kA/gTY'^L^/^ [0, we get the above resuh 
for Lc from the condition 5x{Lc) ~ W . 

The full shape of the PDFs is determined by the rel- 
ative cumulants Rp^x = [XP]c/[XY. For X ^ N , x, c 
or xc, our results yield Rp,x ^ X ^ numerical 

coefficient. Therefore, the system parameters enter the 
PDF shapes only through i?2,x, which is interestingly 
independent of the observable X , and given by 



R 



2,X 



(VKL)- 



gT ) WL 



(11) 



The length scales in the final expression are = ag, the 



1/2 3/2 



the mean 



separation of FLs, and = (gT/kA) 
longitudinal distance between collisions of FLs. To ob- 
tain a size independent PDF, anisotropy requires looking 
at finite-size samples of width W ^ L'', with roughness 
exponent ( — 2/3. The resulting PDF is however, still 
non-universal, depending on i?2,x ^ [sT/kAY^'^. 



The above predictions for the scaling of PDFs, or cu- 
mulants, can be tested experimentally by measurements 
on different realizations of randomness, drawn from the 
same distribution of impurities. Generating many such 
different realizations could in fact be quite easy, depend- 
ing on the system under study. For example, in the case 
of the FL lattice of Ref . 0| experiments can be performed 
on the same sample, with different realizations of ran- 
domness generated by simply rotating the sample with 
respect to the external magnetic field. Each (finite-size) 
realization of randomness yields a characteristic value 
for thermodynamic observables, providing a reproducible 
'fingerprint' of the configuration. The statistics of these 
measurements is then described by the cumulants calcu- 
lated here. 
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